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ABSTRACT 

Earth surface motions on the order of one to several cm per year have been ob- 
served. These motions ccHild be due to land subsidance, i.e. the gulf coast of Texas 
subsides about 5 to 10 cm per year; likewise, similar values hold for portions of the 
Florida coast), crustal uplift (i.e. dilatency vdiich is a phenomena observed to 
precede earthquakes primarily along thrust faults), and loading (i.e. due to 
large dams, etc.). Knowledge of these motions is of practical importance to 
government and local agencies. 

A technique is discussed here for the accurate (i.e. to within fractions of cm/yr) 
detection of these motions utilizing the latest space technology, ft is shown that 
over a six day period and assuming a 50% cloud cover (i.e. as experienced over 
the last few years of laser operation) utilizing spaceborne precision ranging 
systems, intersite distances on the order of 5 to 15 km (dependent mostly on the 
beam width of the laser) can be determined in the vertical and horizontal compo- 
nents with errors in the 0.5 to 1.5 cm range. These errors are almost independent 
of ground survey errors up to 0.25 meters and orbit errors up to 200 meters. 

A spaceborne laser ranging system is assumed to range simultaneously to two or 
more ground emplaced retroreflectors. The fundamental advantage derived from 
simultaneous ranging is the elimination to first order of errors due to the system. 

This means elimination of (a) bias errors in the ranging system, (b) errors due to 
propagation effects, and (c) errors associated with the spacecraft's motion in its orbit. 

In conclusion, it can be stated that horizontal and vertical errors in relative 
distance determination up to 25 km are expected to be in the 0.5 to 1.5 cm range 
using only a 6 day mission. Ihis assumes very modest errors in the ranging 
noise cr of 5 cm (present GSFC laser system precision) , orbital errors of 200 
meters, and a priori knowledge in the intersite distances of 0.25 meters. 

These results make this system very attractive as a new tool for monitoring 
very small earth surface motions. 
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INTRODUCTION 

Accurate determination of intersite distances has been pursued at Goddard for 
many years. This work began in 1964 when the first satellite (BE-B) equipped 
with laser corner reflectors was launched. Continuous improvement of ground 
based laser ranging systems made it possible to determine intersite distances 
of a thousand km to within perhaps tens of centimeters. As early as 1968, J. 
Rosenberg suggested reversing the system; that is, using corner cubes on the 
ground and precision laser systems in orbit. Serious consideration was given 
to such a system during 1971/72 at Goddard. At this time work began to design 
a precision laser ranging system for the space shuttle laboratory. Note that 
the following analysis applies to both electronic transponders as well as laser 
corner cube retroflectors. The analysis considers only ranges which are in- 
dependent of the systems used for their determination. An analysis using pos- 
sible rai^e rates will be published at a later date. 

One impetus came from the possible use of the application of such a system 
for earthquake prediction. C. Scholz, L. Sykes and Y. Aggarwal presented a paper at 
the AGU meeting in Washington, D. C. (Ref. 1) in the spring of 1973 stating that dila- 
tancy, among other phenomena, is a realprecursor to earthquakes. As outlined in 
their paper, the ground rises around an active zone. A well documented case is, to 
quote the example used by Scholz, the 1964 Niigata earthquake (M = 7.5)* (See 
Figure 1). Measurements Indicated rather clearly a constant vertical motion of 
the ground of about 12 cm over 60 years (0.2 cm/yr). This is followed by a quite 
rapid change of 5 cm over 3 years (1.7 cm/yr). Subsequently, the motion stops 
for 3 to 4 years till the earthquake occurred. A system as described here could 
perform such measurements. Accurate (cm-level) relative intersite distances 
are needed in many other practical applications. Such a system could further be 
installed to monitor small motions near or on large dams, major construction 
sites, shore facilities and structures, etc. 


I. SYSTEMS CONCEPT 

As was shown by Scholz, et al. (1073) Ref. 1, the region of rock dilatancy extends 
one fault length (L) on each side of an active fault zone (total length = 3L) and one 

*M s Eorthquake magn'tude as measured on Richter Scale. 
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Figure 1. (Reference 1) 
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length on each side (width = 2L). Thus, if one wishes to observe the rise of the 
ground, an array of ground transponder cubes which covers this active area 
must be used. (A = 6L^) as shown schematically in Figure 2. 

Ihe smallness of the motions to be observed over several years requires the 
use of a rather specific ranging and/or tracking systems concept. In essence, 
one wants to determine very accurately (cm-range) a distance 1 D ( or better, its 
changes between two ground points using satellite technology. In general, there 
are three major obstacles or error sources to overcome if one wishes to com- 
pute extremely accurate vector distances of 5 to 20 km, that is accurate to 0.5 
to 1.5 cm from a spacecraft. These are: 

1. Orbital errors of the spacecraft 

2. Bias errors in the ranging system 

3. Atmospheric propagation errors 

These can be eliminated to first order by using range differencing. Practically, 
this means sending one pulse from the spacecraft to two or more ground stations, 
subtracting these measurements and using only their "differences” to compute 
I D[ , the distance between neighboring ground stations (i.e. transponders or laser 
comer cubes). It is thus evident by studying Figure 3 that the above mentioned 
error sources will not play any major role (only second order) since 


- (Pfej + (^kj+l 

- (Pkj “ ^kj +1 ^ 


For bias errors in the ranging system 


LASER (RADIO) 
BEACON SPACECRAFT 


REFLECTORS 
TRANSPONDER) 
IN BEAM 



Figure 2. Rock Dilatancy Determination from Space 
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because both the and (j + 1 )‘^ reflectors* are interrogated by the same 
signal (within the beam). For orbital errors and atmospheric refraction, 


+ l ^ (second order terms) 

In summary, using range differences as the basic "measured" quantity elimi- 
nates all first order errors in the determination of the distance | 6 |. 

Using one pulse to cover more than one corner cube (or transponder) unfortunately 
raises a power problem. From a signal to noise ratio vantage point it is preferable 
to send a single beam out at a time to one ground station. This does, however, 
create a problem by introducing the effects of errors 1 and 2 mentioned 
before. For instance, a two second time interval between pulses means their 
origin (spacecraft position) has separated in space by say 15 km. This means 
the orbit error introduced by those 15 km along the orbit pass will now increase 
the error in the determination of j d[. The bias error in range, from firing to 
one cube and the subsequent firing to the other cube will further directly in- 
fluence the determination of [ D| . Thus the control of bias errors in the ranging 
system and the orbit error becomes an important systems design factor. A 
detailed analysis of this second possibility will subsequently be published. 


II. MATHEMATICAL MODELS 

1. DIRECT BASELINE ESTIMATION USING SIMULTANEOUS RANGING 

During the estimation of the accuracy of the intersite vector distance ! D j from 
scalar range difference measurements, (p^^ - + l ), two coordinate systems 

are employed. The first is a local topocentric system, Z 3 x 1 , centered at the 
observer. The second is a geocentric earth fixed system, Uj^ j , to w^ich space- 
craft motion is related. The intersite vector distance !d| is defined in the 
same Z 3 X 1 system in which a surveyor would work. The variation in |d| is the 
quantity needed to determine small (cm-range) horizontal and vertical dis- 
placements in the earth’s upper crust. 


If ground transponders ore used, this is not so simple, since the deloy time in each tronsponder 
is not the same. (1 nanosec = 30 cm). 
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As defined in Ref. 2 (Kshn, Vonbun,1966), the Z -coordinate system is centered 
at the dt)server with the 

- axis directed towards local east 
z 2 -axis directed towards local north 
Zj -axis directed along normal to local horizon plane 

In order to relate the local ^ coordinate system with the dynamics of a 
spacecraft, a transformation into a geocentric earth fixed system must be made, 
lids system, named the U-coordinate system, is defined as follows: 

U| -axis directed towards Greenwich Meridian 

Uj -axis normal to (u^ , U 3 ) axes 

U 3 -axis directed along earth's axis of rotation. 

TTie transformation needed later from the Z ( 3 ^ ^ j -coordinate system into tiie 
U-coordinate system is as follows (see reference 2): 


I 


Z(3k1) = Ri(^/2-0)R3(77/2 +\) [11(3^ 


*(3x 1)J 


(1) 


Mtdmre 


X : is the position vector of the observer In the U-system 
o I ) 


Rj(7t/2 -cp) ^ 
R3(7t/2 +\) • 


and 

are rotation matrices. 


The vector 0 is now to be computed from rai^e difference measurements which 
are scalar quantities. It is to be shown how D is computed. The basic measure- 
ment, slant range p , from the ground statiems to the spacecraft or vice versa 
can be expressed in terms of the coordinate systems defined in equation 1 in 
matrix notation as follows: 


p = (ZTZ)‘/2 


( 2 ) 


where 


and 


Z^Z = U^U + S^S - ZS^u 

= geocentric distance to spacecraft 
(S^S)’'"^ = geocentric distance to ground station 
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By assuming that at time t^ , a signal Is sent out by the ranging system onboard 
the spacecraft aiKl received at the and (j 1)*^ station. The range at each 
station Is by (2) expressed as follows: 


P., = UU^U). t (S^S). - 2SJUJ’''’ (3a) 

p.,.1 = ItU'U), . (S^S),„ - 2sy„uj’'’ (3b) 


with similar expressions as (3a) aixl (3b) resulting for times 

^k+l * ^+2‘ * * *’ 

The range difference measurement is obtained by subtracting equation 3b from 
3a which then represents the fundamental observation equation, that is; 




- [(U^U), + (STS).„ - 


(4) 


In arriving at equation 4, the assumption made is that both stations, the j aixi 
(j + 1)‘^ simultaneously observe the laser pulse at time t^. 

Since the range difference measurement is subject to errors as all measure- 
ments are, and since D Is derived from these measurements, D Is also subject 
to errors. We estimate the errors of D due to errors in the range difference 
measurements (p^. )* To do so, standard linear estimation theory 

(Ref. 3) is u&sd for a set of k observations of the same cubes. That is, by vary- 
ing (4) and using first order terms of the Taylor expansion 


* 1)" ^(kx6)*0 ’*^®(kK3)*j + ^'(k* 3) *j + l 

(6x1) (3* n ( 3 k 


one obtains: 


+ *(kxl) 


1 ) 
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v^re 


y(kxl)“ ^ ^kj " ^ki ♦ l^^k* i; 


(kx6) 


j± [iv-s,r- -i_ 

iTkj ^kj + 1 J 


A. 


/ 

3Uo 


Xo ^ 5Uo 
(6x 1 ) (6K 1) 


B, , * - — [U^ - s.]^ 

<kx3) p k J 
'^kj 


s. . SS, 
(3» 1) 


Vm ■*«,.. 

(3x 1) 

€ s vector of observation errors. 


2. INTERSITE DISTANCE DETERMINATION AND ITS ACCURACY 

An estimation of the intersite distance between two ground emplaced comer 
cube and their errors is obtained from range diffemnce measurements using 
the standard least squares estimiUion technique (Ref. 3, 6. 7). In order to 
minimize the effects due to geopotential uni^rtainties, multiple short orbital 
arcs of 1 to 1-1/2 revolutions are used. In the least squares process, e.stimates 
are to be c^tained simultaneously for each orbital arc's state vector as well as 
the coordinates of the comer cul^s. The intersite distances and their errors 
are computed by a least squares solutton. A mathematical descrij^ion of the 
estimation process now follovra. 

The generalized matrix equation which represents all the observation equations 
obtained from simultaneous ranging from the spacecraft to two ground emplaced 
laser retroreflectors reads as follows: 
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( 6 ) 


1 ‘ {6r+6)**^(«r*6)«! + | 


where 



r * number of different orbital arcs for which the state vector 
is to be estimated r = 1, 2, . . . p; 

P 2 1 


^ ^ kj ; kj 

» • 1 


total number of observations associated 
with each orbital arc. 


i > > > (6r + 6) 


t(kjK 1) 



<k^Kl) 


■C X 1 


Vector of all measurement residuals 
observed during the entire *SPEAR* 
data collection period. 


r. 


'i* (6r *6) 


*1 

t 

: • 

« 

H, 








IW|, ,) 


• ’ 


( 
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*1 

‘ j 
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(k^, j ^ 







‘ j > 


• 

• 

• 

4, 

. ». 

: 'S.i' 

-j 

(6r ♦ 6) 


Matrix of sensi> 
tivity coefficients 
associated with each 
measurement. 
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X 


1 

0 

(6x 1) 


T 

*(6r +6)x 1 


( 6 >‘ 1 ) 


S 


i (3x 1) 


(3x 1) 


(6 r +6)x 1 


Vector of error sources which are 
major contributors to the measui'ement 
residuals. These error sources are 
the states associated with each orbital 
arc used for calculation of measured 
minus calculated range differences 
(i.e. residuals) and corrections to the 
comer cube coordinates. This vector 
is estimated by the least square process. 


(k,xl) 




X 1 


(k^xx). 


(^xl) 


Vector of residuals which reflect the 
neglect of second order terms of Taylor 
expansion made for representing the 
errors in the measurements. 


The estimate of the vector T is obtained from -i range difference measurements 
witii > > 1 and a priori information about the state of each orbital arc and 
station survey. That is; 


T(6,.6)x, = r^*'T + w;n ' [nw-'Y+»;jT„] 

(6f+6)x(6r + 6) (6r+6)xl 


( 7 ) 


vdiere 

W"* - A priori covariance matrix associated with system accuracy. 

= A priori covariance matrix associated with orbit accuracy and 
° station survey for each retroreflector. 
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» 


The errors associated with the estimate of ^ are obtained from the covari- 
ance matrix given by the following equations 


E(TTT) 


ErTw-ir 



-1 

(6r + 6)x (6r +6) 


(8) 


Equations (7) and (8) are now used to (a) compute the intei::.ite distance and (b) 
its errors. The intersite distance, that is, the difference between the comer 
cube position vectors (s. - Sj+i ) is obtained by a coordinate transformation 
of the vector! . This transformation reads as follows; 


D (Sj “ Sj + i)3x 1 - Qsx (6r + 6)^(6r + 6)*l 


where 

[®(3x6r)|l(3x3)l “ ^( 3 x 3 )] 3 ^ ^ ^ ^ ^ g ^ 

0 = a n\jl I matr ix 

1 = the identity matrix 

The error in the intersite distance is now given by the covariance matrix of 
the intersite distances, that is: 


E[(», - 3, = (1») 


which upon introducing Equation 8 finally reads: 


E[(s, -^h)’' 1 ( 3 . 3 , = (0[r"w-'r (11) 
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Equation 11 representing the intersite distance error, is the basic equation used 
for the SPEAR System error analysis. Numerical results using (11) are depicted 
in Figures 4 & 5. 

m. ERROR ANALYSIS RESULTS 

In order to evaluate the accuracy with which intersite distances can be deter- 
mined using the SFEAR concept, a parametric error analysis study was per- 
formed utilizing the mathematical model described by Equation 11. For this 
study very moderate systems errors were assumed. The results obtained, are 
presented in Figures 4 and 5. All assumptions concerning the source of errors 
and their magnitude are indicated on each figure. 

Although orbital errors up to 200 meters were used in the^error analysis, their 
contribution to the horizontal and vertical components of D consistently re- 
mained negligible. Furthermore, survey errors associated with the location of 
the retroreflectors also do not significantly change the errors in the vertical 
and horizontal components of D . This means that no extra care has to be taken 
in the emplacement of ground equipment. 

Figure 5 shows the dependency of the vertical and horizontal distance errors on 
the noise of the ranging systems used. As can be seen, a 5 cm system is ade- 
quate, producing errors in the 0.7 and 1.3 cm ranges. As stated earlier, 5 cm 
ground laser ranging systems are now operational at Goddard Space Flight 
Center, thus no difficulties are foreseen in achieving this precision for a space- 
borne laser system. 


IV. PRACTICAL APPLICATIONS 

The SPEAR system can be developed as a payload for the shuttle applications 
program. As a matter of fact, a breadboard model is under development at 
Goddard (Ref. 5). 

Such a system can easily be utilized in detecting small relative variations of 
the Earth’s upper crust in the 0.3 to 0.5 cm/yr range. Monitoring, is thus 
possible, of land subsidance as it occurs at coastal regions of Florida and Texas, 
as well as construction sites such as dams and even large buildings. Dilatency, 
the vertical uplift along seismic zones, can further be followed and thus studied 
as a precursory effect which takes place before many shallow foci earthquakes 
(Ref. 1, 3, 4). The San Andreas fault area and the Niigata area in Japan are 
examples. 
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INTERSITE DISTANCE ERROR (cm) 


\ LASER NOISE; 5cm 

I ORBIT ERRORS: UP TO 200 METERS 

I \ SURVEY ERRORS: UP TO 0.2 METERS 

1 \ ZENITH ANGLE: 70® 

\ \ SAMPLE RATE: 1 MEAS/SEC 



BOTH ERROR COMPONENTS ARE ALMOST 
LINEAR FUNCTIONS OF THE SYSTEM 
NOISE ONLY. (DOUBLING LASER SYSTEM 
NOISE DOUBLES THESE ERRORS) 
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Figure 4. Intersite Distance Errors vs. Time in Orbit 
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ORBIT ERROR: 200 METERS 
SURVEY ERROR: 0.2 METERS 
BASELINE SEPARATION: 25km 
ORBIT ALT.: 400km 
ORBIT INCL.: 50° 

DAYS IN ORBIT: 6 
ZENITH ANGLE: 70° 

SAMPLE RATE: 1 MEAS/SEC 


LASER SYSTEM UNCERTAINTY (cm) 


Figure 5. Intersite Distance Errors vs. Laser System Uncertainty 
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V. CONCLUSIONS 


This paper has shown that the SPEAR system can be used to determine inter- 
site distances up to 15 km depending only on the beam width or beam splitting 
of a spaceborne laser system within a 6 day shuttle mission to a precision 0.5 cm 
to 1.5 cm (assuming 50% cloud coverage). As anticipated, (a) survey errors (up 
to 0.25 m), (b) orbital errors (up to 200 m), and (c) range bias errors (many 
meters) play only a minor role, making the system potentially useful for many 
practical applications where very small relative motions are to be monitored. 


REFERENCES 

1. Scholz, C. H., Sykes, L. R., and Aggrawal, Y. P. "Earthquake Prediction: 

A Physical Basis"; Science, Vol. 181, Number 4102, August 1973. 

2. Kahn, W. D., Vonbun, F. O., "Tracking Systems, Their Mathematical Models 
and Their Errors" Part II: Least Squares Treatment, NASA TN D-3776; 
December 1966. 

3. Deutsch, Ralph; "Estimation Theory" Prentice Hall, 1965. 

4. Siry, J. W., "Crustal Motion Measurement by Means of Satellite Techniques", 
X-590-73-273; May 1973. 

5. Minott, P. O., "Design of a Ground Based Laser Retroreflector Array for 
Space to Ground Geodesy", Unpublished Report, Revised March 1975. 

6. Old Dominion Systems, Inc., "User's Guide to Data Preparation, Navigation 
Analysis Program (NAP 3.1)," April 1974. 

7. J. J. Lynn, "NAPCOV Program Documentation". Old Dominion Systems, 

Inc. November 1973 (revised May 1, 1974). 


16 


